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Abstract 

The electron-phonon interaction and related transport properties are investigated in monolayer 
silicene and M0S2 by using a density functional theory calculation combined with a full-band 
Monte Carlo analysis. In the case of silicene, the results illustrate that the out-of-plane acoustic 
phonon mode may play the dominant role unlike its close relative — graphene. The small energy of 
this phonon mode, originating from the weak sp 2 it bonding between Si atoms, contributes to the 
high scattering rate and significant degradation in electron transport. In M0S2, the longitudinal 
acoustic phonons show the strongest interaction with electrons. The key factor in this material 
appears to be the Q valleys located between the T and K points in the first Brillouin zone as 
they introduce additional intervalley scattering. The analysis also reveals the potential impact of 
extrinsic screening by other carriers and/or adjacent materials. Subsequent decrease in the actual 
scattering rate can be drastic, warranting careful consideration. Finally, the effective deformation 
potential constants are extracted for all relevant intrinsic electron-phonon scattering processes in 
both materials. 

PACS numbers: 71.15.Mb, 72.10.Di, 72.20.Ht, 73.50.Dn 
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I. INTRODUCTION 



Very recently, the attention to low- dimensional materials has expanded beyond the best 
known example of this kind — grapheneP^ In particular, silicene^] and molybdenum disul- 
fide^^ have gained much interest due to their unique properties in electronics, optoelectron- 
ics and magnetics. Silicene is expected to share certain superior properties of graphene due 
to its structural similarity and the close position in the periodic table. More importantly, 
it is compatible with the current silicon-based technology and can be grown on a number 
of different substrates On the other hand, atomically thin M0S2 is a semiconductor 
with a finite band gap that ranges from approximately 1.3 eV to 1.9 eV depending on the 
thicknessP It has been used as the channel material in field effect transistors with promising 
resultsP^E^ In addition, monolayer M0S2 offers the possibilities of interesting spin and valley 
physics utilizing the strong spin-orbit couplingP^^ 

Characterization of electronic transport, particularly the intrinsic properties, is critical 
for assessing and understanding the potential significance of a material. In the case of sil- 
icene, many of the crucial parameters are presently unknown due to the brief history of 
this material. In comparison, notable advances have been made in M0S2 lately. Exper- 
imental investigation of transistor characteristics deduced the channel mobilities ranging 
from ~200 cm 2 /Vs to ~1000 cm 2 /Vs at room temperature,^^ while a theoretical study 
estimated an intrinsic phonon- limited value of ~410 cm 2 /Vs based on a first-principles cal- 
culation of electron-phonon interaction.^ However, questions remain regarding the intrinsic 
electron transport in M0S2. For instance, those extracted from transistor current-voltage 
(I-V) measurements are indirect accounts and can be strongly affected by extrinsic factors. 
Similarly, the latter worlP^ includes only the electronic states in the conduction band min- 
ima at the K points in the first Brillouin zone (FBZ); the impact of Q valleys located along 
the T-K symmetry directions (which correspond to the minima of bulk M0S2) were not 
considered. A detailed investigation is clearly called for. 

In this paper, we examine intrinsic transport properties of monolayer silicene and M0S2 
by taking advantage of first-principles analysis and full-band Monte Carlo simulation. Along 
with the electronic band structure, the phonon spectra and electron-phonon coupling ma- 
trix elements are calculated for all phonon branches within the density functional theory 
(DFT) formalism.* 21 * 22 ! The obtained electron scattering rates are subsequently used in the 
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Boltzmann transport equation to compute the intrinsic velocity-field characteristics with a 
full-band Monte Carlo treatment. The calculation results are compared with the available 
data in the literature and the key factors affecting electron transport in these materials 
elucidated. The investigation also provides the effective deformation potential constants 
extracted from the first-principles results. 

II. THEORETICAL MODEL 

Both monolayer silicene and M0S2 are hexagonal crystals. To account for their delicate 
atomic structures accurately, the calculations are performed in the DFT framework, as it 
is implemented in the QUANTUM-ESPRESSO package,^ using ultrasoft pseudopotentials. 
A minimum of 35 Ry is used for the energy cut-off in the plane wave expansion along with 
the charge truncation ~15 times larger. The generalized gradient approximation is used 
for the exchange-correlation potential for silicene, while the local density approximation is 
adopted for M0S2. The momentum space is sampled on a 36x36x1 Monkhorst-Pack grid 
with no offset (silicene) or on a 18x18x1 grid (M0S2) for electronic band calculation. The 
simulated cells are optimized in both cases until the atomic forces decrease to values less 
than 0.015 eV/A. 

Each phonon is treated as a perturbation of the self-consistent potential, created by all 
electrons and ions of the system, within the linear response [i.e., density functional pertur- 
bation theory (DFPT)]. 21 The calculation of the potential change due to this perturbation 
gives the value of the electron-phonon interaction matrix element^ 



where \i, k) is the Bloch electron eigenstate with the wavevector k, band index i, and energy 



atomic displacement associated with the phonon from branch v with the wavevector q and 
frequency oo v ^; and M is the atomic mass. Numerical calculations of lattice dynamics are 
conducted on a 18x18x1 Monkhorst-Pack grid. Indices i, j are dropped hereinafter as only 
the first (lowest) conduction band is considered. 

With the electron-phonon interaction matrix from the first-principles calculation, the 




(1) 



Eiy, AV^ SCF is the derivative of the self-consistent Kohn-Sham potential with respect to 



3 



scattering rate of an electron at state |k) can be obtained by using Fermi's golden rule, 

- = IT 52 IgqkH^q^k+q " ^,q " #k) + (JV„,q + l)5(^k-q + ^*,q " #k)], (2) 
q,« 

where iV V) q = [exp^c^.q/^B^) + l]" 1 is the phonon occupation number, ks the Boltzmann 
constant, and T the temperature. As we are interested in the intrinsic scattering probability 
that is not limited to a specific carrier distribution (and thus, the Fermi level), our formula- 
tion assumes that all final electronic states are available for transition (i.e., nondegenerate) 
in the bands under consideration. 

For transport properties, a Monte Carlo approach with full-band treatment is adopted. 
All of the details described above, including the wavevector (k,q) dependence of the scatter- 
ing matrix elements [e.g., Eq. Q], are accounted for. This allows solution of the Boltzmann 
transport equation beyond the conventional relaxation time approximation. 

III. RESULTS AND DISCUSSION 
A. Monolayer Silicene 

Earlier first-principles studies have shown that the stable structure for monolayer silicene 
has a low-buckled configuration.^ While planar and high-buckled cases lead to imaginary 
phonon frequencies around the V point indicating an unstable structure, the low-buckled 
construction provides well separated phonon branches and positive frequencies. The origin 
of buckled geometry in silicene is the weakened ir bonding of the electrons in the outer 
shell. Compared with graphene, which has very strong ix bonding and planar geometry, 
the interatomic bonding distance is much larger in silicene, which decreases the overlap of 
p z orbitals and dehybridizes the sp 2 states. Accordingly, the planar geometry cannot be 
maintained. In our analysis, the lattice constant a is optimized to be 3.87 A, with the 
buckling distance of 0.44 A, in good agreement with Ref. HI 

Figure [T] shows the outcome of electronic and phononic band calculation in monolayer 
silicene. The Fermi velocity extracted from the Dirac cone is around 5.8 x 10 7 cm/s that is 
roughly one half of that in graphene [see bands 7r and n* in Fig. [jja)]. While this result is 
in agreement with other theoretical predictions j 24 l 25 l a value as high as 1 x 10 8 cm/s was also 
claimed in the literature P As for phonons in Fig. [lib), six branches are identified with two 
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atoms per unit cell. The transverse acoustic (TA) and longitudinal acoustic (LA) phonon 
dispersion relations are well approximated by sound velocities in the long- wavelength limit; 
t>TA = 5-4 x 10 5 cm/s and v L a = 8-8 x 10 5 cm/s. Although the out-of-plane acoustic (ZA) 
phonon exhibits an approximate q 2 dependence near the center of the Brillouin zone, its 
sound velocity can also be estimated; vzk = 6.3 x 10 4 cm/s. An interesting point to note 
in the phonon dispersion is the discontinuities in the frequency derivative of the highest 
optical branch that, similar to graphene, appear at the high symmetry points, T and K. 
These discontinuities are referred to as Kohn anomalies] 24 * 26 ^ induced by unusual screening 
of lattice vibrations by conduction electrons. Sharp cusps typically indicate strong electron- 
phonon coupling. Calculated phonon energies at the T, M, and K points in the FBZ are 
summarized in Table [H 

The electron-phonon interaction matrix elements obtained for the electron at the Dirac 
point [i.e., k = (47r/3a, 0)] are plotted in Fig. [2] as a function of phonon wavevector q. Kohn 
anomalies, illustrated by the three peaks at three equivalent K points in the transverse 
optical (TO) mode and another at the zone center for the longitudinal optical (LO) branch, 
are not as distinct as those observed in graphene.^ Overall, coupling of optical phonons with 
electrons appears to be relatively weak. In comparison, the acoustic phonons show much 
stronger interaction. Particularly striking is the large strength of ZA phonon coupling, 
unlike in graphene. Due to the buckled geometry (originating from the weak it bonding 
mentioned earlier), silicene does not maintain certain key characteristics of ideal planar 
lattice, especially the reflection symmetry with respect to the atomic plane. As such, the 
symmetry consideration, in which only the in-plane phonons can couple linearly to two- 
dimensional (2D) electrons,^ no longer applies. An increased role of ZA phonons is clearly 
expected. 

The scattering rates calculated at room temperature (T = 300 K) are shown in Fig. [3j 
The result is plotted specifically for electrons with wavevector k along the K~r direction. 
Since the integration in Eq. ^ is over the entire FBZ, both intravalley (K — > K) and 
intervalley (K — > K') transition events are included. As the interaction matrix elements 
illustrated above suggest, acoustic phonons have much larger scattering rates than optical 
modes. Specifically, the ZA branch provides the dominant contribution, which can be at- 
tributed to the observed large coupling strength as well as the small phonon energy near the 
zone center (i.e., a large occupation number AzA,q)- This also indicates that the scattering 
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rates are very sensitive to the phonon energies (or equivalently the value of uza)- Since an 
accurate description of ZA dispersion in the long wavelength limit requires a well converged 
calculation with a sufficiently dense grid, care must be taken when evaluating accuracy of 
the data in the literatureP^ 

Figure H(a) provides the drift velocity versus electric field at different temperatures ob- 
tained by full-band Monte Carlo simulations. The intrinsic mobility estimated from the 
figure is approximately 1200 cm 2 /Vs and the saturation velocity (defined at 100 kV/cm) 
3.9 x 10 6 cm/s at 300 K. When the temperature decreases, both the mobility and the sat- 
uration velocity enhance due to the suppression of phonon excitation; the respective values 
at 50 K are 3.0 xlO 4 cm 2 /Vs and 6.2 x 10 6 cm/s. The drift velocities show a slight negative 
slope at high fields that becomes more pronounced at low temperatures. This phenomenon 
(i.e., the negative differential resistance) can be explained, at least in part, by the nonlinear 
band dispersion at high electron energies as in grapheneP^ 

The calculation results discussed above demonstrate the intrinsic properties of silicene. 
When this material is synthesized or placed on a substrate, however, additional scattering 
sources such as surface polar phonons and impurities must be considered, which could de- 
grade the performance further. A topic that may need additional attention is the role of 
ZA phonons in the presence of a supporting material. As recent measurement of graphene 
in-plane thermal conductivity attests,^ even a weak binding between a 2D crystal and the 
substrate could dampen the ZA vibrations substantially. Moreover, it is reasonable to an- 
ticipate that the extent of this suppression would be dependent on the detailed interaction 
between two materials. Since ZA phonons provide the dominant role in the electron-phonon 
interaction in silicene, it (i.e., the damped vibration) could actually lead to sizable reduc- 
tion in the scattering rate. To gauge the impact, transport characteristics are also studied 
without the ZA scattering. As shown in Fig. |4](b), the mobility experiences an increase of 
greater than threefold (3900 cm/s), while the saturation velocity goes up more modestly 
(5.6 x 10 6 cm/s). This estimate may be considered an upper limit for silicene on a substrate. 

B. Monolayer MoS 2 

In the present DFT calculation for monolayer M0S2, the optimized lattice constant is 
3.13 A, consistent with other theoretical studies.^ 1 30 1 31 1 Furthermore, this value is in good 
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agreement with 3.15 A determined experimentally in bulk MoS2.^ The resulting electronic 
and phononic band dispersion is depicted in Fig. |5j As shown, monolayer M0S2 is a semi- 
conductor with a direct gap of 1.86 eV at the K point — a number within a few percent 
from a recent measurement of 1.9 eVP Our calculation also predicts the presence of second 
energy minima only about 70 meV higher. These so-called Q valleys are located along the 
T-K axes at approximately the half-way points [e.g., Q=(0.34 x 27r/a,0) vs. K=(47r/3a, 0)]. 
At present, the energy separation -Eqk between the K and Q valleys is unsettled with the 
estimates ranging from 50 meV to 200 meVP 311 Since this is a crucial parameter for electron 
transport, a more extended discussion is provided later in the paper in relation to intervalley 
scattering. The band dispersion relations around the energy minima are nearly quadratic 
and can be well described by the effective mass approximation. For the K valleys (i.e., 
the energy minima at the equivalent K points), the extracted longitudinal and transverse 
effective masses are almost identical, m l K = m\r = 0.50mo- On the other hand, the Q valleys 
yield mq = 0.62mo, itlq = 1.0m with the longitudinal direction defined along the T-K axis. 
m denotes the electron rest mass. 

Monolayer M0S2 has the symmetry of the point group D 3 h, with nine branches of phonons. 
The irreducible representations associated with each phonon mode, together with the polar- 
ization (longitudinal or transverse), help denote all of the vibrational modesj^EESI as plotted 
in Fig. |5j The E" modes are degenerate at the T point. These two modes are the in-plane 
optical vibrations, with two S atoms moving out of phase and Mo atom static. The E' 
modes are polar LO and TO phonons, with Mo atom and two S atoms moving out of phase. 
Due to the coupling with the macroscopic electric field, there is LO-TO separation at the F 
point, which slightly lifts the LO(E') mode upward on energy scale. A non-analytical part 
is added to the dynamic matrix resulting in a small splitting of about 0.3 meV [not visible 
due to the energy scale of Fig. |5^b)]. The Ai and A 2 ' branches are two out-of-plane optical 
phonon modes. Ai is also referred to as the homopolar mode, with two S atoms moving 
out of phase and Mo atom static. In the Ag mode, Mo atom and two S atoms vibrate out 
of phase. The three lowest branches are LA, TA, and ZA modes, with sound velocities of 
^la = 6.6 x 10 5 cm/s, t>TA = 4.3 x 10 5 cm/s. The phonon energies at different symmetric 
points are summarized in Table [TTJ 

Figures [6] and [7] show the electron-phonon interaction matrix elements for the initial 
electron state at k = K [=(47r/3a, 0)] and k = Q [« (2ir/3a, 0)], respectively, for TA, 



7 



LA, TO(E'), LO(E'), and Ai (or homopolar) phonon modes. The contribution from the 
remaining four branches is found to be negligible due to the weak coupling. The matrix 
elements for k = K demonstrate a threefold rotational symmetry (i.e., 120°), while those of 
k = Q show the reflection symmetry with respect to the q x axis. As expected, the LO(E') 
phonons near the T point possess the characteristics of Frohlich coupling through the induced 
macroscopic electric field typical of polar materials. Since the relative electronic potential 
is not periodic in the long wavelength limit DFPT does not yield a correct value to the 
electron-phonon interaction matrix. For an approximation, the matrix element of LO(E') is 
interpolated at F by using the values from the nearby q points. This [i.e., LO(E')] and A x are 
the only two modes that have non-zero scattering matrix as q — > (intravalley scattering); in 
the other three branches, the matrix elements only have first-order components, |g k | ~ q, 
leading to |g q ^. ,kl ~~ ^ 0- With regard to intervalley scattering that requires large q phonons, 
a number of different transition processes are possible as shown in Figs. [6]T) and[7](f). For 
instance, Fig. [6] indicates strong electron-phonon interaction at the symmetry points M in 
the phonon momentum space (denoted as q = M for simplicity) for all modes except LO(E'); 
these phonons can induce electron transition from K to Q' valleys. Another example is the 
phonons at q = K 7 for all five modes in Fig. [7J which can be associated with electron 
scattering from Qi to Q4. 

The electron-phonon scattering rates are calculated as a function of electron energy using 
Fermi's golden rule. Figure [8] gives the rates for electrons in the K valleys at room temper- 
ature, while the result for Q-valley electrons is shown in Fig. [9] Similarly to silicene, the 
wavevector k of the initial electronic state is chosen along the K-r or Q~r axis, respectively. 
As can be seen from the figures, the LA mode provides the largest scattering rates consistent 
with its large coupling strength. The discontinuities or steps in the curves represent either 
the onset of optical phonon emission or intervalley scattering. For instance, the abrupt in- 
crease observed in the rate of LA phonons at ~100 meV in Fig. [8]^a) can be attributed to 
the above mentioned strong K — > Q' transition via emission of a LA phonon with q = M. 
Since this phonon energy is approximately 30 meV (see Fig. [5J, the final state energy of 
~70 meV indeed matches to the Q-K separation £qk- Similarly, the onset of transition via 
absorption can be found around 40 meV in Fig. (8|b). The large density of states in the Q 
valleys (evident from the large effective masses) makes the contribution of this scattering 
even more prominent. If, on the other hand, all of the final states in the Q valleys are 
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excluded, the total scattering rate for the K-valley electron reduces drastically to approx- 
imately 2 x 10 13 s -1 , which is consistent with the prediction of an earlier first-principles 
calculation.^ The observed difference of an order of magnitude clearly illustrates the strong 
dependence of the scattering rates on Eqk- The inconsistency of this value in the recent 
publications^^ adds difficulty to accurately evaluating the role of Q valleys. 

Utilizing the scattering rates, the velocity vs. field relation is obtained by a full-band 



Monte Carlo simulation at four different temperatures. As shown in Fig. [TOj the extracted 
mobility decreases from 4000 cm 2 /Vs at 50 K to about 130 cm 2 /Vs at room temperature 
while the saturation velocity changes from 7.6 x 10 6 cm/s to 3.4 x 10 6 cm/s. The small 
mobility and saturation velocity can be attributed to strong electron-phonon scattering as 
well as the heavy effective masses. With massive electrons that hinder acceleration and 
many states to scatter into (e.g., K, K', Q, Q' valleys), this is an expected outcome. 

Compared to a recent theoretical estimation^ of 410 cm 2 /Vs and experimentally ex- 
tracted values^ as high as 1090 cm 2 /Vs, however, our mobility is significantly smaller, 
requiring a careful analysis of the discrepancy. Two factors are identified that could yield at 
least a partial explanation. First, let us examine the issues surrounding the Q-K separation. 
With inconsistencies reported in several first-principles results on this sensitive quantity 
(see the discussion above), it is not unreasonable to imagine that our DFT calculation may 
have also experienced similar inaccuracies. If -Eqk proves to be substantially larger than 
the estimated 70 meV, then the Q valleys would have a limited influence on the low-field 
mobility and can be neglected in the calculation as in Ref. [201 (with 200 meV). In this case, 
our simulation estimates the K-valley dominated mobility of 320 cm 2 /Vs that is essentially 
in agreement with the earlier prediction (410 cm 2 /Vs)P*l Clearly, both first-principles mod- 
els produce a consistent picture of K-valley electron dynamics including intrinsic scattering 
with relevant phonon modes. The difference is the relative significance of Q valleys (e.g., 
70 meV vs. 200 meV). As such, further clarification of intrinsic mobility in monolayer M0S2 
may need to be preceded by accurate experimental determination of Eqk- 

Even when the influence of Q valleys becomes negligible, a sizable disparity remains be- 
tween the theoretically obtained mobility and the highest value claimed experimentally^ 
(e.g., 320—410 cm 2 /Vs vs. 1090 cm 2 /Vs). This is puzzling since the theoretical estimates 
of the intrinsic parameters are generally expected to provide an upper limit to the mea- 
sured data that tend to experience additional mobility degrading, due to external scattering 
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sources. When examining Ref. [15j however, it becomes clear that the extracted data is 
not a direct measurement of intrinsic mobility; rather, it is the channel mobility that is 
strongly influenced by the details of the structure and the bias conditions. A particularly 
salient point is that the deduction relied on the transistor I-V data measured when the M0S2 
channel is populated by electrons (i.e., degenerate). This condition deviates from the under- 
lying assumption of nondegeneracy in our calculation, necessitating consideration of, among 
others, the screening effect. In a low-dimensional system, it is known that the screening 
due to the degenerate electrons can lower the bare scattering rates substantially (thus, af- 
fecting the mobility) j 34 l 35 l as previously demonstrated experimentally in the AlGaAs/GaAs 
structures.^ The screening effect can be included by renormalizing the electron-phonon 
interaction through the dielectric function; i.e., gq k — > gq ik /e(g). | 34 | a5 1 

To gauge the potential significance in monolayer M0S2, we adopt a simple model for the 
dielectric function based on Thomas- Fermi screening of only the K- valley electrons: i.e., 
e(q) = 1 + Qtf/q, where the screening wavevector q TF = Am^e 2 /U 2 k. Here, the factor 
of 4 accounts for the spin and valley degeneracies, e is the electron charge, and k is the 
background dielectric constant. Subsequent calculation with a rough estimate of e(q) shows 
that the scattering rates can experience a decrease of about an order of magnitude through 
screening. A corresponding increase of the mobility is estimated to be well over 1000 cm 2 /Vs 
that is more consistent with the value extracted from the transistor I-V characteristics.^ 
Thus, it is evident that the screening must be taken into account accurately when the carrier 
density becomes degenerate in M0S2. A detailed analysis of e(q) is, however, beyond the 
scope of this investigation as our focus is on the properties of intrinsic electron-phonon 
interaction. 

C. Deformation Potential Model 

For practical applications, it would be convenient to approximate the ab initio results for 
electron-phonon coupling by a simple analytical model. Particularly useful in the present 
case is the deformation potential approximation. Under this treatment, the coupling matrix 
(j, k + q| A V^ SCF |2, k) in Eq. (JT| can be expressed in the first order (D^q), or in the zeroth 
order (Dq)^ The first-order deformation potential constant {D\) is adopted to represent 
the coupling matrices for the acoustic phonon modes in the long wavelength limit (i.e., 



10 



intravalley scattering). In comparison, those involving the near zone-edge acoustic phonons 
(i.e., intervalley scattering) are treated by using the zero f,l -order deformation potential (-Do) 
in a manner analogous to the optical modes. In the latter case (Do), the phonon energy is 
assumed independent of the momentum for simplicity. The obtained analytical expressions 
of the scattering rates are then matched to the first-principles results by fitting the effective 
deformation potential constants. 

For silicene, the intravalley scattering rate by acoustic mode v (= LA, TA, ZA) is obtained 

as 

D\k B T 



1 

1(1) 



£ k . (3) 



Here p is the mass density (= 7.2 x 10~ 8 g/cm 2 ) and v v denotes the sound velocity, for 
which we can take the value of v%x = 6.3 x 10 4 cm/s, wta = 5.4 x 10 5 cm/s, and wla = 
8.8 x 10 5 cm/s, respectively, as discussed earlier. On the other hand, the rate of optical 
phonon scattering (both intravalley and intervalley transitions) as well as the intervalley 
acoustic phonon scattering can be expressed by the following form: 



'k,« 



g 2h 2 Vppu v 



D 2 

" [(E k + hLo v )N v + (E k -hcu v )(N v + l)e(E k -hco v )}, (4) 



where 9(x) is the Heavyside step function and N v = [exp(hu v /kBT) + for phonon 
mode v. Since the phonon dispersion is treated constant in Eq. (|4]), the intravalley optical 
phonon scattering approximates the zone-center (i.e., T) phonon energies for hu v (v = LO, 
TO, ZO). In the case of intervalley scattering via acoustic or optical phonons, hu v takes 
the respective phonon energy at the zone-edge K point corresponding to electron transition 
K -H- K'. The specific values used in the calculation can be found in Table [Tj When matched 
to the first-principles rates, the effective deformation potential constants can be extracted 
for each scattering process as summarized in Table |III| A particularly interesting point to 
note from the result is that all three acoustic phonons show generally comparable values of 
D\ and D despite the large scattering rate of ZA mode (see also Fig. [3]). Clearly, this mode 
(ZA) couples strongly with electrons but not enough to prevail over other acoustic branches; 
its dominant contribution is due mainly to the small phonon energy as discussed earlier. 

The electron-phonon scattering processes in M0S2 are much more complicated as the 
deformation potentials need to be determined for both K- and Q-valley electrons. While 
feasible, it is not practically useful to define the effective interaction constants individually 
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based on the mode of involved phonons and the transition types. Accordingly, we adopt 
a simplified description by combining the appropriate contributions into just two modes, 
acoustic and optical, respectively. 

Using the effective mass approximation for the band structure near the valley minima, 
the scattering rate for intravalley acoustic phonon scattering (i.e., K — > K or Qi — > Qi by 
both LA and TA phonons; see Figs. [6] and [7|) is given by 



mtD(k B T 



Mos 2 s 



(5) 



r (i) 

where p = 3.1 x 1(T 7 g/cm 2 for MoS 2 and m* is the density-of-states effective mass for the 



K or Q valley (final state), m* = ymj,mj, (p = K,Q). For electrons in the Q valleys, strictly 
speaking, an additional factor of 0(i?k — -Eqk) is multiplied to the left side of Eq. ^ to 
account for the non-zero energy at the bottom of the Q valley. By taking the sound velocity 
v s = vla (= 6.6 x 10 5 cm/s), the value of D\ is estimated to be 4.5 eV and 2.8 eV in the K 
and Q valleys, respectively. 

The analytical expression that describes intravalley and intervalley optical phonon scat- 
tering as well as intervalley acoustic phonon scattering rate is obtained as 



r (2) 



m* D 2 

= dd-^iN^ + (N v + 1)A 2 ], (6) 

MoS 2 



where g& is the valley degeneracy for the final electron states, and Ai and A 2 denote the on- 
set of scattering for phonon absorption and emission, respectively. For instance, Ai = 1 and 
A 2 = Q(Ek — hu) v ) for electron transitions between the K valleys, whereas Ai = O^k - -Eqk) 
and A 2 = Q(Ei c — hu v — -Eqk) for transitions between the Q valleys. The factors correspond- 
ing to intervalley transfer between K and Q valleys can be constructed accordingly, where 
-Eqk may be treated as a potentially adjustable parameter. Tables IV and [V] summarize 



the initial/final electron states, the phonon momentum that is involved (in the form of its 
location in the momentum space), and the extracted deformation potential constants for 
each transition process considered in the investigation. For a given phonon momentum, the 
actual value hu v used in the analytical calculation is the average of the relevant phonon 
modes. Specifically, the acoustic (optical) phonon energy is obtained as the average of LA 
and TA [TO(E'), LO(E'), and Ax] modes at the respective symmetry points given in Table pi} 
An additional point to note is that the estimate of Dq P at the T point includes the effect 
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of Frohlich scattering by the LO(E') mode. 20 While this is a mechanism physically distinct 
from the deformation potential interaction and must be handled separately, its impact is rel- 
atively modest, at least for the electrons in the K valley. Accordingly, the present treatment 
is considered adequate. Further simplification of the model may also be possible judging 
from the narrow range of values in Dq c and Dq P (mostly in the low to mid 10 8 eV/cm). 

IV. SUMMARY 

We have performed a first-principles calculation together with a full-band Monte Carlo 
analysis to examine electron-phonon interaction and the intrinsic transport properties in 
monolayer silicene and M0S2. The results clearly elucidate the role of different branches 
as well as the intr a/inter- valley scattering. The predicted intrinsic mobility for silicene is 
approximately 1200 cm 2 /Vs, with saturation velocity of 3.9 x 10 6 cm/s at room temperature. 
In the case of M0S2, the K-valley dominated mobility gives approximately 320 cm 2 /Vs, while 
the intrinsic value reduces to about 130 cm 2 /Vs when the energy separation of 70 meV is 
used between the K and Q minima. The estimated saturation velocity is 3.7 x 10 6 cm/s. The 
investigation also illustrates the significance of extrinsic screening, particularly in numerical 
evaluation of transport characteristics. The extracted deformation potential constants may 
prove to be useful in further studies of these materials. 
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TABLE I: Phonon energies (in units of meV) at the symmetry points for monolayer silicene. 



Phonon modes 


r 


K 


M 


ZA 





13.2 


13.0 


TA 





23.7 


13.4 


LA 





13.2 


13.5 


ZO 


22.7 


50.6 


50.7 


TO 


68.8 


50.6 


56.7 


LO 


68.8 


61.7 


64.4 
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TABLE II: Phonon energies (in units of meV) for TA, LA, TO(E'), LO(E'), and Ai (or homopolar) 
modes at the T, K, M and Q points in the FBZ of monolayer M0S2. 



Phonon modes 


r 


K 


M 


Q 


TA 





23.1 


19.2 


17.9 


LA 





29.1 


29.2 


23.6 


TO(E') 


48.6 


46.4 


48.2 


48.0 


LO(E') 


48.9 


42.2 


44.3 


44.2 


Ai 


50.9 


51.9 


50.1 


52.2 
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TABLE III: Extracted deformation potential constants for electron-phonon interaction in silicene. 



Phonon mode 


Intravalley 


Intervalley 


ZA 


2.0 eV 


6.1 xlO 7 eV/cm 


TA 


8.7 eV 


1.4xl0 8 eV/cm 


LA 


3.2 eV 


4.2xl0 7 eV/cm 


ZO 


6.3xl0 7 eV/cm 


4.3 xlO 7 eV/cm 


TO 


1.8xl0 8 eV/cm 


1.4xl0 8 eV/cm 


LO 


1.9xl0 8 eV/cm 


1.7xl0 8 eV/cm 
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TABLE IV: Extracted deformation potential constants for electron-phonon interaction in M0S2 for 
electrons in the K valley [see also Fig. [6^f)]. 



Phonon 


Elect 


ron 




Deformation 


Momentum 


Transition 




potentials 


r 


K -> 


K 


n° p - 


Df =4.5 eV 
=5.8 x 10 8 eV/cm 


K' 


K -> 


K' 


T) ac - 


=1.4 x 10 8 eV/cm 
=2.0 x 10 8 eV/cm 


Q' 


K -> 


Q 


T) ac - 
iy Q - 


=9.3 x 10 7 eV/cm 
= 1.9 x 10 8 eV/cm 


M 


K -»• 


Q' 


T) ac - 

n op - 


=4.4 x 10 8 eV/cm 
=5.6 x 10 8 eV/cm 
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TABLE V: Extracted deformation potential constants for electron-phonon interaction in M0S2 for 
electrons in the Qi valley [see also Fig. [71(f)]. Multiple equivalent valleys for the final state specify 
the degeneracy factor larger than one in Eq. 



Phonon 


Electron 




Deformation 


Momentum 


Transition 




potentials 


r 


Qi 


-+Qi 




Df =2.8 eV 
=7.1 x 10 8 eV/cm 


Q 3 (Q 5 ) 


Qi -> 


Q 2 (Q 6 ) 


T) ac - 

n op 


=2.1 x 10 8 eV/cm 
=4.8 x 10 8 eV/cm 


M 3 (M 4 ) 


Qi -> 


Q 3 (Q 5 ) 


T) ac - 

n op - 


=2.0 x 10 8 eV/cm 
=4.0 x 10 8 eV/cm 


K' 


Qi 


-»• Q4 


T) ac - 
D op 


=4.8 x 10 s eV/cm 
=6.5 x 10 8 eV/cm 


Qi 


Qi 


-»• K 


T) ac - 


=1.5 x 10 8 eV/cm 
=2.4 x 10 8 eV/cm 


M 2 (M 5 ) 


Qi 


-►K' 


T) ac - 


=1.5 x 10 8 eV/cm 
=2.4 x 10 8 eV/cm 



21 




22 



(a) ZA 



(b) ZO 




-0.5 0.5 
(c) TA 



-0.5 0.5 
(d) TO 



0.5 



e 



-0.5 



0.5 



CM 



-0.5 



-0.5 0.5 
(e) LA 



-0.5 0.5 
(!) LO 



0.20 



0.15 



0.10 



0.05 



0.0 



-0.5 0.5 
cjx (2jt/a) 



-0.5 0.5 
qx(2rc/a) 



FIG. 2: (Color online) Electron-phonon interaction matrix elements |gk +q jj (in units of eV) from 
the DFPT calculation in silicene for k at the conduction-band minimum K point [i.e., (47r/3a, 0)] 
as a function of phonon wavevector q for all six modes v. 
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FIG. 3: (Color online) Electron scattering rates in silicene via (a) emission and (b) absorption of 
phonons calculated at room temperature. The electron wavevector k is assumed to be along the 
K-r axis. 
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FIG. 4: (Color online) Drift velocity versus electric field in monolayer silicene obtained from a 
Monte Carlo simulation at different temperatures: 50 K (square), 100 K (triangle), 200 K (dia- 
mond), and 300 K (circle). The results in (a) consider the scattering by ZA phonons, while those 
in (b) do not. 
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FIG. 5: Electronic and phononic band structures of monolayer M0S2 along the symmetry directions 
in the FBZ. The conduction-band minimum at the K point serves as the reference of energy scale 
for electrons. 
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(a) TA (b) LA 




-0.5 0.5 -0.5 0.5 




FIG. 6: (Color online) (a)-(e) Electron-phonon interaction matrix elements |gk +q iJ (in units of 
eV) from the DFPT calculation in M0S2 for k at the conduction-band minimum K point [i.e., 
(47r/3a, 0)] as a function of phonon wavevector q. Only the branches with significant contribution 
are plotted; i.e., TA, LA, TO(E'), LO(E'), and Ai (or homopolar) modes, (f) Schematic illustration 
of intervalley scattering for electrons in the K valley. 
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FIG. 7: (Color online) (a)-(e) Electron-phonon interaction matrix elements |gj^. + k | (in units of 
eV) from the DFPT calculation in M0S2 for k at the Q point [i.e., Qi ~ (27r/3a, 0)] as a function 
of phonon wavevector q. Only the branches with significant contribution are plotted; i.e., TA, LA, 
TO(E'), LO(E'), and Ai (or homopolar) modes, (f) Schematic illustration of intervalley scattering 
for electrons in the Q valleys. 
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FIG. 8: (Color online) Scattering rates of K-valley electrons in M0S2 via (a) emission and (b) 
absorption of phonons calculated at room temperature. The electron wavevector k is assumed to 
be along the K-r axis. 
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TO(P)_ _ . . : 
(b) 

0.10 0.15 0.20 0.10 0.15 0.20 

Electron energy (eV) 

FIG. 9: (Color online) Scattering rates of Q-valley electrons in M0S2 via (a) emission and (b) 
absorption of phonons calculated at room temperature. The Q-K separation energy I?qk (= 70 
meV) denotes the onset of curves as the K- valley minimum serves as the reference (zero) of energy 
scale. The electron wavevector k is assumed to be along the G-r axis. 
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FIG. 10: (Color online) Drift velocity versus electric field in monolayer M0S2 obtained from a 
Monte Carlo simulation at different temperatures with -Eqk = 70 meV. When electron transfer to 
the Q valleys is not considered, the mobility increases to approx. 320 cm 2 /Vs at 300 K. 
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